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We investigate the kinetics of A + _B ^ reaction with the local attractive interaction between 
opposite species in one spatial dimension. The attractive interaction leads to isotropic diffusions 
inside segregated single species domains, and accelerates the reactions of opposite species at the 
domain boundaries. At equal initial densities of A and B, we analytically and numerically show 
that the density of particles (p), the size of domains {€), the distance between the closest neighbor 
of same species {Iaa), and the distance between adjacent opposite species {£ab) scale in time as 

1/3 1/3 2/3 

p ^ t , £aa ~ t , and £ ~ £ab ~ t ' respectively. These dynamical exponents form a new 
universality class distinguished from the class of uniformly driven systems of hard-core particles. 

PACS numbers: 05.40.+j, 82.20.-w 



I The irreversible two-species annihilation reaction A + 
B ~^ has been intensively and widely investigated as 
a basic model of various phenomena in physics P, Q], 
chemistry Q , and biology Q ■ The reaction includes two 
species A and B which are initially distributed at ran- 

] dom in space. The reaction of opposite species instan- 
taneously takes place with a rate k when two particles 
of opposite species encounter on the same site (generally 
within a reaction radius) during the motion of particles. 
The reaction forms the third inert species, which is then 
disregarded in the motion of A and B. 

' For the same initial density of A and B, pa{0) — Pb{Q), 
the mean-field equations predict that pA and pB decay 
linearly in time as {kt)~^. However it turned out that the 
random fluctuation of the initial number of two species 
evolves to make a segregation into A— or i?— rich area 

! Qi IE IE Q' The fluctuation and segregation develop 
in time, so that the reactions of two opposite species 

' take place only at the boundaries of two adjacent seg- 
regated domains. As a result, in sufhciently low dimen- 
sions, the effect of fluctuation leads to the anomalous 
kinetics. That is, the evolution of the density of particles 
strongly depends on fluctuations, and cannot be derived 
from mean-field rate equations. 

The density decay has been known to depend on the 
motion of the particles and the mutual statistics of par- 
ticles. For isotropic diffusions, the density pA{t) scales 

[ in time t as PAit) ~ t"''/* in d-dimensions (d < dc = 4) 
IE IE, II II ll 03 ■ The dimension 4 is the upper criti- 

[ cal dimension above which the density decay follows the 
mean-field rate equations, so p{t) ^ for d > dc- With 
the global relative drift of one species, pA{t) scales as 
PA{t) - for d < 3 0. The hard-core (HC) in- 

• teraction between identical particles is irrelevant to the 
case of the isotropic diffusion and the relative drift 0. 
However when both species are uniformly driven to the 
same direction, the HC interaction completely changes 
the asymptotic scaling as pA ^ t^^'^+^y^ for d < 2, t"'^/'^ 
for 2 < d < 4, and finally t^^ for d > 4 [H E^. With- 
out the HC interaction, pA{t) decays as pA ~ t"''/* as in 
the isotropic diffusion due to Galilean invariance. Recent 
studies on the reaction under Levy mixing 111 and on 



scale free networks 14j showed that some mixing mech- 
anisms that homogenizes reactants can suppress the role 
of the fluctuations. 

In reality in which oppositely charged particles re- 
combinate into inert particles such as electron-hole re- 
combinations in photoluminescence ^El and particle- 
antiparticle reactions in the early universe 0, the local 
attractive interaction between opposite species should be 
much more important than the global uniform biases. In 
this paper, we investigate the kinetics of A + £? ^ re- 
action with the attractive interaction between opposite 
species in one dimension, which may model composite 
systems of oppositely charged particles much more effec- 
tively than the model with the global and uniform drift. 
When a particle is surrounded by two same species neigh- 
bors such as BAB, the central particle (^4) performs ran- 
dom walks. If two opposite species particles surround a 
particle such as AAB, the central particle (A) is ballisti- 
cally driven to its opposite species {B). As a result, the 
attractive interaction depends on the local configurations 
of adjacent particles, and accelerate the reactions of op- 
posite species at the boundaries of segregated domains. 
However inside segregated domains, each particle has the 
same neighboring species, so the motion is isotropic diffu- 
sion. Hence the situation belongs to neither the relative 
drift nor the isotropic diffusion. Due to the isotropic 
diffusion inside domains, the HC interaction should be 
irrelevant in this case. 

With the local attractive interaction in one spatial di- 
mension, we analytically and numerically show, regard- 
less of the existence of the HC interaction, that the 
density of particles (p), the distance between the clos- 
est neighbor of same species (Iaa), the size of domains 
(£), and the distance between adjacent opposite species 
{£ab) scale in time as p ^ t~^/^, Iaa ^ t^^^, and 
(■ ~ f-AB ^ t^^^ respectively. These dynamical expo- 
nents form a new universality class distinguished from 
the class of uniformly driven systems of hard-core parti- 
cles 0, 01 , where the exponents for Iaa , ^ and £ab are 
still controversial. The two features of ballistic and diffu- 
sive motions result in pentagonal space-time trajectories 
of bulk particles (Fig. 1(a)), which allow us to derive the 
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asymptotic scaling analytically. 

We consider a configuration in which A and B species 
are randomly distributed on an one dimensional lat- 
tice with an equal initial density, pa(0) = p_b(0). A 
randomly-chosen particle performs either isotropic or bi- 
ased random walks depending on the configurations of 
neighboring particles. When the chosen particle is sur- 
rounded by two same species neighbors such as BAB, 
the chosen particle (A) performs isotropic random walks. 
If two opposite species particles surround a particle such 
as AAB, the chosen particle (A) is constantly driven to 
the its opposite species (B). 

In the region of a length £, the number of A species 
is initially — /9^(0)^ ± \/pA(fi)i and the same for 
Nb- After a time t ^ particles travel throughout 
the whole of the region, and annihilate in pairs. The 
residual particle number is the number fluctuation in the 
region so we have the relation V( or pA ~ 

for a given length ^ H, 0| ■ the processes evolve, the 
system becomes a homogeneous collection of alternating 
A-rich and i?-rich domains. To characterize the structure 
of segregated domains, we introduce three length scales 
as in Ref. j^. The length of the domain (i) is defined 
as the distance between the first particles of adjacent 
opposite species domains 0. The length Iab is defined 
as the distance between two adjacent particles of opposite 
species, while ^aa(^bs) is the distance between adjacent 
A{B) particles in a A{B) domain. These length scales 
asymptotically increase in time as 

With the attractive interaction, as a bulk particle inside 
single species domains diffuses isotropically until it be- 
comes a boundary particle, its space-time trajectory is 
pentagonal as shown in Fig. 1(a). These pentagonal tra- 
jectories should be self-similar (self-affine) fractal struc- 
tures, because they should have the scaling symmetry 
under the scaling transformation x' = bx (space-domain 
scaling) and t' — b^t (time-domain scaling) with 6 > 1 
due to the power-law behavior in Eq. (1). A typical base 
unit of the self-similar pentagonal trajectories of adjacent 
opposite domains are schematically depicted in Fig. 1(b). 
This base unit allows us to calculate a time Tg needed to 
remove the unit of the space-domain size £ surrounded 
by one scale larger ones. Then the size of the larger unit 
increases by £ during Tg so we have 

d£/dt - £/Te , (2) 

which gives the dynamic exponent z. 

As only boundary particles of each domain have two 
opposite species neighbors, the boundary particles bal- 
listically annihilate due to the attractive interaction. It 
takes a time ri — £ab for two boundary particles to an- 
nihilate in the base unit. The second particle from the 
boundary isotropically diffuses during time ti until the 
boundary particles annihilate. After the time ti, the sec- 
ond particle becomes a boundary particle, and constantly 




(a) (b) 

FIG. 1: (a) A snapshot of space-time trajectories of y4 -f 
B with the local attractive interaction between opposite 
species, (b) The magnified schematic space-time trajecto- 
ries of one pentagonal base configuration of adjacent opposite 
species domains. Subscripts {1,2,. ..,n} indicate the order of 
the positions of particles from a given domain boundary. 

moves to its counterpart during t2 = £ab + £aa- So it 
takes time T2 = ti + t2 in total for the second particle 
to annihilate. Similarly, the nth particle from the ini- 
tial boundary will annihilate after Tn = Tn-i + tn, where 
tn = ^AB + [n — 1)£aa for n > 2. From the recurrence 
relation of r„, we find 

Tn = n£AB + n{n - \)£aaI2 . (3) 

As the number of particles in a domain of size £ is A'^ ~ 
V^, the time needed to annihilate the domain in the 
base unit is given by 

n ~ N,£ab + iV/ tAA , (4) 

for Ng ^ \. In above calculations, we consider the 
mean positions of bulk particles, and neglect the in- 
crease of £AA(fy by diffusions during the annihilation of 
the base unit. After a smaller unit is completely annihi- 
lated by being embodied into larger one, the remainder 
of particles redistribute over the larger unit increased by 
the size of the annihilated unit. Hence we approximate 
lAAit) = ■■■ = £AA{t + T„) = • • • = £AA(t + Tf) during 
the annihilation of the unit. 

The scaling of £aa is simply £aa from the rela- 

tion lAAit) ~ 1/ p{t). Hence £aa scales as £aa ~ t^^"^^ 
with ZAA = 2z. On the other hand, the change of £ab 
during T£ should be order of Ng £aa because of Ng succes- 
sive annihilations of two opposite particles at boundaries. 
So we get MAs/dt ^ A£ab /t^ ^ £/t. With the scaHng of 
t ^ £^, we find that £ab follows the same scaling as £, i.e. 
^AB ^ t^^^ with ZAB — z. Finally using the relations, 
£aa ^ V£, Ng ^ ^/l, and £ab ^ ^, we find ti from Eq. 
(4) 

n - e''' ■ (5) 

Substituting Eq. (5) into Eq. (2) and integrating the 
resultant equation, we finally arrive the following scaling 
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relation 



(6) 



As a result, we find z = zab =3/2 with the attractive 
interaction between opposite species. 

From the scaling of £, £ab, and Iaa, asymptotic decays 
of various densities can be extracted. The density of total 
particles {p = 2pa), the density of pairs of adjacent same 
species {paa — Pbb), and the density of pairs of adjacent 
opposite species {pab) scale as 



t PAA ~ t~ 



PAB ~ t" 



(7) 



As p is order of p ^ 1/%/^, we have p ^ t ^l'^'^ with 
a — l/2z = 1/3. PAA follows the same scaling of 
p due to Paa ^ ^I^aa ~ 1/V^ so paa ~ with 
aAA — a — 1/3. Finally pab is Pab ~ 1/^, which leads 
to PAB ^ t^^l^ with aAB = ^/z = 2/3. Using self- 
similar structures of space-time trajectories and scaling 
arguments for fluctuations of 0, IE > following 
exponents for the reactions A + B ^ with the local 
attractive interaction between opposite species 



a = aAA = 1/3 , aAB = 2/3 , 

Z = ZAB = 3/2 , ZAA = 3 . 



(8) 



Intriguingly and incidentally p{t) decays with the same 
exponent 1/3 as that of the uniformly driven hard-core 
particles [H in which the driven motion of a sin- 
gle species domain was arg ued to be described by the 
noisy Burgers equation 'll]. However for the attractive 
interaction case, the | decay of p{t) comes from the in- 
terplay of isotropic diffusions inside domains and ballistic 
annihilations at boundaries. For scalings of inter-domain 
distances and others, our results of (8) are completely dif- 
ferent from those of [l2|, where i ~ t^^^^,iAB t^^^, and 
^AA ~ t^^^- Hence we conclude that the local attractive 
interaction between opposite species forms a new univer- 
sality class of irreversible A + B ^ reactions. 

To see the validity of our analytic results, we now 
want to report the simulation results for the model with 
the attractive interactions. With equal initial density of 
PAiO) — Pb(0), a and B particles distribute randomly on 
a lattice of size L. In the simulations we consider both 
hard-core (HC) particles and the particles without the 
HC interactions, which we call the bosonic particles. In 
the model with HC particles there can be at most one 
particle of a given species on a site. In the bosonic model 
there can be many particles of the same species on a site. 
As we shall see, the simulation results are independent 
of the HC interactions. 

All the simulations are done on the one-dimensional 
chains with the size up to L = 3 x 10^ and the initial 
densities are always taken as pa{0) = Pb{Q) = 0.2. We 
average p(t), pAA{t), and pAsit) up to 10^ time steps over 
7200 independent runs. In Fig. 2, we plot the densities 
and their effective exponents defined as 
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FIG. 2: densities (inset) and effective exponents for the densi- 
ties in the model with hard-core particles (a) and those in the 
model with bosonic particle (b). From top to bottom, each 
line corresponds to oaa, a and aAB- In the inset, the order 
is p, pAA, and pab respectively. 
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FIG. 3: The scaling collapse plot of measured p{t) for Eq. 
(10). Used sizes of lattices are L = 2^*, 2^* and used par- 
ticles are hard-core particles. We set a = 0.333 and z — 1.5. 
The inset shows the raw data for p{t). We do not show the 
scaling collapse plot for bosonic particles, which are almost 
identical to the main plot. 



and similarly for others. As you can see in Fig. 2, 
the data for HC particles (Fig. 2(a)) are almost identi- 
cal to those for bosonic particles (Fig. 2(b)). While aAB 
still shows larger fluctuations for both HC particles and 
bosonic particles, a and aAA nicely converge to the same 
value. We estimate a ~ 0.33(1), aAA — 0.33(1), and 
ctAB = 0.68(2) for both HC and bosonic models, which 
agree well with the predictions (8). 

To estimate the dynamic exponent z, we measure the 
densities for various L from 2^** up to 2^*. With the 
scaling assumption 



p{L,t)^t-^f{t/L^) 



(10) 



a(t)=ln[p(t)/p(</2)]/ln2. 



(9) 



and the estimate a = 0.333, we observe the best data 
collapse at z = 1.50(2) which also agree well with the 
prediction of (8) as shown in Fig. 3. 

For the time dependence of average distances defined 
in Eq. (1), we measure £, (-aa, and £ab under the same 
measurement conditions as those of densities. The effec- 
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FIG. 4: The average distances(inset) and effective exponents 
for the distances in the model with HC particles (a) and with 
bosonic particles (b). From top to bottom, each line corre- 
sponds to 1/z, 1/zab, and 1/zaa respectively. In the inset 
the order of lines is the same as in the main plot. 

tive exponents of the distances are defined similarly as in 
Eq. (9) except the minus sign. The results for the average 
distances are shown in Fig. 4. Here we can also see that 
the data for HC particles (Fig. 4(a)) are almost identi- 
cal to those for bosonic particles (Fig. 4(b)). From Fig. 
4, we estimate 1/z = 0.683(3), 1/zaa = 0.338(3), and 
^/zAB = 0.66(1) respectively for both HC and bosonic 
particles. 

All the simulation results in Figs. 2-4 numerically con- 
firms the predictions (8). They also confirms the irrele- 
vance of the HC interactions in the case of the attractive 
interaction unlike that of the uniformly-driven systems 
Is, 11]. The isotropic diffusion inside segregated domains 
leads to Galilean invariance of domains so the asymp- 
totic scaling behavior is not affected by the existence of 
hard-core interaction. 

In summary we investigate the anomalous kinetics of 
A + B ^ reactions with the attractive interaction be- 



tween opposite species in one dimension. As reactions 
proceed, initial random fluctuations in the particle num- 
ber develop. As a result. A- and i3-rich domains appears 
alternatively, and annihilations of opposite species take 
place only at the boundaries of the closest neighboring 
domains as in ordinary A + B ^ reactions 0, 0, 0] ■ 
However the reactions at domain boundaries are acceler- 
ated by the attractive interaction, while particles inside 
domains isotropically diffuses until they become bound- 
ary particles. The interplay of isotropic diffusions of bulk 
particles and ballistic annihilations of boundary particles 
lead to pentagonal self-similar trajectories which allow 
us to analytically calculate the dynamic exponent z and 
others as in Eq. (8). We numerically confirm the pre- 
dictions (8) by means of Monte Carlo simulations. The 
results are irrelevant to the existence of the hard-core 
interaction between same species particles. 

The anomalous density decay of appear to be- 

long to the same universality class as uniformly driven 
systems of hard-core particles, which was argued to be 
described by noisy Burgers equations Tl^. However our 
system is not in the same universality class as the uni- 
formly driven system, because scaling behaviors of basic 
distances are different from each other 111. The dif- 
ference can infer from the underlying mechanisms. In 
our model, isotropic diffusions inside domains lead to 
Galilean invariance of domains so the hard-core inter- 
action has no effects on the reactions. Furthermore there 
is no global bias to one direction which change the ki- 
netics of hard-core particles. Only boundary particles 
feel bias to opposite species, which is the essential phys- 
ical factor and distinguishes our model from the models 
in Refs. P, 0, 0- We conclude that the attractive in- 
teraction between opposite species is the key feature of 
the new universality class characterized by exponents in 
Eq. (8), and it is another path to the anomalous density 
decay of t^^/^. 
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